FIG. 1. Figurel f^(q 2 ) with s = 33GeV 2 ,m b = A.lGeV 
FIG. 2. Figure2 /^(g 2 ) in LO WF with a s corrections 
FIG. 3. Figure3 f^(q 2 ) in LO WF with q s =0 
FIG. 4. Figure4 f^ w (q 2 ) with s = 33Ge V 2 ,m b = 4.7GeV; a s corrections. 
FIG. 5. Figure5 f^(q 2 ) with a2(jjt) = 0.12, a4(fi) = 
FIG. 6. Figure6 fg (0) as function of Borel Parameter M 2 
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Figure2 f + B7r (q ) in LO WF with a s corrections 
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Figure3 1\ n (q ) in LO WF with a=0 
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B - 7T WEAK FORM FACTOR WITH CHIRAL CURRENT IN THE LIGHT-CONE SUM RULES 
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In this article, we calculate the B — » it transition form factor fg^{q 2 ) by including perturbative 
0(q s ) corrections to the twist-2 terms with chiral current in the light-cone QCD sum rule approach. 
The corrections to the product /s/^(g 2 ) in the leading twist approximation is found to be about 
30%, while a similar magnitude corresponding to 0(a s ) corrections for /b(? 2 ) in the two-point sum 
rule cancel them and result in small net corrections for fg^iq 2 )- Our results confirm the observations 
made in previous light-cone QCD sum rule studies. 

CO ■ PACS numbers 13.20.He ll.55.Hx 
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I. INTRODUCTION 

Quantum chromodynamics (QCD) is the appropriate theory for describing the strong interaction at high energy 
region, however, the strong gauge coupling at low energy destroys the perturbative expansion method. The long 
distance properties of QCD, especially the hadronic matrix elements can provide many important information for 



m 
>. 

, understanding and testing the standard model and beyond. The exclusive semileptonic decay B — > irlvi can be used 

to determine the CKM parameter | V^t, | [1]. However, it requires a reliable calculation of the form factor fg^{q 2 ) 

defined by (n(p) | b-f^u | B(p + q)) = 2f^ 7T (q 2 )p l , + (/^(g 2 ) + /^(f? 2 ))^, with p and p + q being the tt- and B-meson 

^v^j [ four-momentum, respectively. fg^(q 2 ) plays a negligible role for semileptonic decays into the light leptons I = e,/i. 

In Ref. [2], the authors propose a formula called QCD factorization approach for B — > 7T7t, ttK and 7rD to deal with 

' nonleptonic decays of B meson. In this approach, the decay amplitudes are expressed in terms of the semileptonic form 
Qh' 

factors, hadronic light-cone distribution functions and hard-scattering amplitudes. The semileptonic form factors, the 

Qh' 

' light-cone distribution functions are taken as input parameters and the hard-scattering amplitudes are calculated by 
' perturbative QCD. Again, the precise knowledge of heavy-to-light form factors plays crucial roles. Among the existing 
approaches, such as QCD sum rules, chiral perturbation theory , heavy quark effective theory and phenomenological 
quark models, the QCD light-cone sum rules (LCSR) approach is very prominent for calculating fg^iq 2 ) [3-5]. 

The light-cone QCD sum rule approach carries out operator product expansion near the light-cone i 2 «0 instead 
of the short distance x ~ while the nonperturbative matrices are parameterized by light-cone wave functions which 
classified according to their twist instead of the vacuum condensates. For detailed discussion of this method , one can 
see Ref. [6]. The LCSR for fg^{q 2 ) is valid at small and intermediate momentum transfer squared q 2 < mg — 2mQX , 
where % is a typical hadronic scale of roughly 500 MeV and independent of the heavy quark mass tuq. 

In this paper, we calculate the form factor fg w (q 2 ) (which is different from Refs. [7-10]) up to twist-4 light- 
cone functions by including perturbative a s -corrections for twist-2 terms using chiral current. Remarkably, the main 
uncertainties of the light-cone sum rules come from the light-cone wave functions. The chiral current approach has a 
striking advantage that the twist-3 light-cone functions which are not known as well as the twist-2 light-cone functions 
eliminated and supposed to provide results with less uncertainties [11]. In fact, only the twist-2 wave function, which 
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is dominant in contributions to the sum rules, has been investigated systematically. The update investigation of 
the twist-3 wave functions can be found in Ref. [12] and the calculations of the form factor fg including the a(s) 
corrections to the twist-3 terms are performed in Ref. [13]. Although the QCD radiative corrections to the twist-2 
term for fg v are proven small in Ref. [8], it is interesting to see whether or not the case for chiral currents. 

The article is organized as follows: correlator and sum rule arc derived in Sec. II; the perturbative correlator are 
calculated to order a s in Sec. Ill; light-cone amplitudes and numerical results are presented in Sec. IV; the section V 
is reserved for conclusion. 

II. CORRELATOR AND SUM RULE 

Let us start with the following definition of B — > 7r weak form factors /s^g 2 ): 

(ir(p)\u^b\B(p + q)) = 2/+ r (g 2 )p, t + [ftM 2 ) + fsM 2 ))^ (1) 

with q being the momentum transfer. Following Ref. [11], we choose a chiral current to calculate the correlator 
function, 

n„(p,?)= i J d 4 xe iqx (n(p)\T{u(x)^(l +>y 5 )b(x), 6(0)*(1 + 75 )d(0)}|0>m 6 , 

= n(rj 2 , (p + qf) Plx + n{q\ (p + qf)q^ (2) 

which is different from that in Ref. [7-10]. Here we take chiral limit p 2 = m 2 = 0. 

We can insert a complete series of intermediate states with the same quantum numbers as the current operator 
bi(l + 75 )d in the correlator to obtain the hadronic representation. After isolating the pole term of the lowest 
pseudoscalar B meson, we get the result: 

H>, q)= Il H (q 2 , (p + qf )p, + U H (q 2 , (p + q) 2 )q^ 

(ir\u^b\B)(B\b l5 d\0)m b v (7r|H7 M (l + 75)|Bg)<Bg|5i(l + 75)d|0)mt 

m 2 B -(p + q) 2 + jf m 2 BH -(p + q) 2 ' U 

The intermediate states Bh contain not only pseudoscalar resonances of the masses greater than tub, but also 
scalar resonances with J p = + , corresponding to the operator bd. Taking into account the definition (B\bij5d\0) = 
mB 2 fB/mb, we obtain: 

Here the contributions of higher resonances and continuum states above the threshold so ar e written in terms of 
dispersion integrations, and the spectral densities p H (s) and p H (s) can be approximated by the quark- hadron duality 
ansatz. We can avoid the pollution from scalar resonances with J p = + by choosing sq near the B meson threshold 
and our final results conform this assumption. 

In the following, we brief out line the calculation of the correlator in QCD theory and work in the large space-like 
momentum regions (p + q) 2 — m 2 < for the bd channel, and q 2 <C m 2 — 0{lGeV 2 ) for the momentum transfer, 
which correspond to the small light-cone distance x 2 w and are required by the validity of the operator product 
expansion method. First, we write down the full 6-quark propagator: 



2 



(0,^)5(0)10)= i J ^e-<*±±^ - J I * 

mt-k 2 



1 k -\- Tfl s~unj / \ 

(5) 



here G^,, is the gluonic field strength, g s denotes the strong coupling constant. Substituting the above b quark 
propagator and the corresponding tt meson light-cone wave functions into Eq.(2) and completing the integrations over 
x and fc, finally we obtain: 



)= 2.Um 2 b f 
Jo 



n( 9 2 , (p + qf)= 2Umi / du 2 n "V . , ,„ + 



<Pw(u) 2ug 2 (u) 



ml — (1 — w)g 2 — u(p + q) 2 (ml — (1 — ?i)g 2 — u(p + q) 2 ) 2 
%ml\g-i{u) + G 2 (u)} f 2< / 9_ L (a i ) + 2<p ± (a>i) - <p\\(oii) - ^(a*) 



(to 2 - (1 - u)q 2 - u(p + q) 2 ) 3 7 1 ml - (1 - ax - ua 3 )q 2 - (ax + ua 3 )(p + q) 2 ) ' 

with C?2(u) = — g 2 (v)dv and .Dc^ = ciaicia2(ia3(5(l — ai — «2 — a 3 ). Here ^ is 7r meson twist-2 light-cone wave 
function, and gi(u), g2{u), ( P±(cti) 7 <p±(ai), <p||(aj), <p|| (ai) are 7r meson twist-4 light-cone wave functions. Their 
detailed expressions are given in section 4. Then we carry out the subtraction procedure of the continuum spectrum 
by the standard procedure and perform the Borcl transformations with respect to (p + q) 2 , and finally obtain the 
result: 



bAQ >~ m 2 B f B e \J A \ u + uM 2 U 3 M 4 



+ 1 dV J Dai 6 (ax tZl) 2 M 2 e ~ mf '" 9 ^ M + ^ ± ^ - ^ (7) 

2 _ 2 

Here A = and sq denotes the subtraction of the continuum from the spectral integral. For technical details, 

one can see Ref. [6,11]. 

III. RADIATIVE CORRECTIONS IN ORDER a s 

In this section, we calculate the perturbative contribution up to a s for twist-2 terms, while the corrections for 
twist-3 terms and beyond are neglected, as they are supposed to be small. Applying Borcl transformation for the a s 
correction terms is tedious , we can facilitate the calculation greatly by writing down the following dispersion integral 
relation: 

I + bM 2 ) = ^ij- f pWtf^e^ds, (8) 

where 

p QCD (q 2 , s) = [ d Uip7T (u)lmT(q 2 , s, u) . (9) 

T JO 

For example, with the zeroth order approximation , one can easily obtain: 

ImT (<2 2 , s, u) = -2nmlS(ml - (1 - u)q 2 - us). (10) 
To order a s , the amplitude can be written as 

T(rx,r 2 ,u)=To(rx,r2,u) + ^^Tx{r 1 ,r 2 ,u) . (11) 
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Here we introduce convenient dimensionless variables r\ — q 2 /m\ and r 2 = (p + q) 2 jm\. There are six Feynman 
diagrams for determining the first order amplitude T\ in pcrturbative expansion, . For simplicity, we perform the 
calculation in Feynman gauge. In calculation, both the ultraviolet and collinear divergences are regularized by 
dimensional rcgularization and renormalized in the MS scheme with totally anticommuting 75. To be more precise, 
the collinear divergences in the hard amplitude are factored out and absorbed in the evolution of the light-cone wave 
function which is determined by the QCD evolution equation [14]. Finally we get the result: 

,2\ 4 



T 1 (r 1 ,r 2 ,u)= 



6(1 + p) 1 , tin, 



1 - In- 



1 



,[(G(p)-G(r 1 )) + (G(p)-G(r 2 ))} 



(l-p) 2 

(^W [G ip) - G{ri)] + T^ [G {P) " ° (r2)] 

2 P+(l-p)ln(l-p) 4 (l- r2 )ln(l-r 2 ) | 2 3-p 



+ 



(1 - u)(n - r 2 ) 



1-P r 2 {l-pY 

(1 - p) In (1 - p) (l-r 2 )ln( l-r 2 )' 
p r 2 



(12) 



with 



f x dt 

p = n +u(r 2 -ri), U 2 (x) = - — ln(l-t), 



G (p) = U 2 (p) + ln'(l - p) - ln(l - p) ( 1 - In ., 

P 



(13) 



As in the calculation of the nondeading order evolution kernel of the wave function ^(u,^), we take the UV renor- 
malization scale and the factorization scale of the collinear divergences to be equal [15-17]. Our results are of the 
same Dirac structure as that of Ref. [8] but with different weight. 

The MS quark mass depends explicitly on the rcnormalization scale p and implicitly on the renormalization 
scheme. A renormalization scheme independent definition of the quark mass within QCD perturbation theory is given 
by the pole mass which is denoted by m£ . As in Ref. [8] , we replace rhb by ml using the well-known one-loop relation: 



m b 



= m£ j 

To 0(a s ), this replacement adds a term, 



1 + 



47T 



,*2 



-4 + 31n- 



4p 



4-31n^ 
P 



(14) 



(15) 



to the renormalized amplitude T\. 

To proceed further according to Eq.(9) we calculate the imaginary part of the hard scattering amplitude for 
r 2 > 1 and r\ < 1: 



!| m/ :„.,,.„„; ±1I^L { 



rn 



*2 



2tt 2 - 6 + 3 In — b - - 2Li 2 (n) + 2Li 2 ( — ) + In' r 2 



r 2 



+2 



(l-r 2 )In(r 2 -l) 01 _ 2 



r 2 



p 2 



ID 



- 2 In (1 -n) +21n(l -n) - 2 ln(l - n) In — b — - 2 In (r 2 - 1) + 21n(r 2 - 1) -21n(r 2 - l)ln— \- 

P 



+o( P -i) 



Hp - 1) 



+ 2 (lnr 2 + 2-21n(r 2 -l) + ln— b - 



r 2 



It' 



P-1 



+ 2- 



1 



r 2 - P \P r 2 
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+2, ( - 21n^l > ) -4^ - 2 ? ^ ( Inp- 2 ln(p - 1) + 1 - In 



{n - r 2 )(r 2 - p) V r 2 ' r 2 -l/ p-1 (ri - r 2 )(p - ri) 
+0(1 -p) 



1 _ m£ 2 \ 1 



2(lnr 2 + --21n(r 2 -l)-ln^|- 

r 2 p 2 / p - 1 



l — n / , ?Ti 



2 1 1-^2 

r 2 - p r 2 



(16) 



lnr 2 + l-21n(r 2 -l)-ln-|- 

(ri - r 2 )(r 2 - p) V P 

here, the operation + is defined by 

F(a;)|+ = lim (f(x)6(1 - x - rj) -8(1 - x - rj) * F(y)dy^ , (17) 

and thus remove the spurious divergences. The above expressions have a little difference compared with the corre- 
sponding ones in Ref. [8] for coefficients of the 6(1 — p) term. 

Substituting Eq.(16) into Eq.(8-9), we can obtain the desired sum rule in 0(a s ) for the form factor in the 
leading twist- 2 approximation. 

IV. LIGHT-CONE AMPLITUDES AND NUMERICAL RESULTS 

Let us choose the input parameters entering the sum rule for fg^(q 2 ) first. To begin with, let us specify the pion 
wave functions. For the leading twist-2 wave function if v (u, p), the asymptotic form is exactly given by perturbative 
QCD <fi v (u,fi — > oo ) = 6u(l — u) [14,18], nonperturbative corrections can be included in a systematic way in term of 
the approximate conformal invariance of QCD and expanded in terms of Gegenbauer polynomials C^ 2 (2u — 1) with 
weight u(l — u). 
To leading order (LO), 

<p„(u, p) = 6u(l - u) £ MPo) ) ° C n /2 ( 2u - !); (18) 

and to non leading order (NLO) [19], 

<p 7r ( U ,p) = 6 U (l-n)^a„(po)ex P (- da 1 -^- J ( C 3 J 2 (2u 1) + ^ £ d k MC 3 k /2 (2u 1) ) , (19) 

with ao = 1. Arguments based on conformal spin expansion allows one to neglect higher terms in this expansion and 
we take n < 4. The coefficients a 2 (po) = 2/3 and 04 (po) = 0.43 at the scale po = 500 MeV have been extracted from 
a two-point QCD sum rule for the moments of <fin(u) [4,18]. The coefficients d^(p) are due to mixing effects, induced 
by the fact that the polynomials C~}J 2 (2u — 1) weight by u(l — u) arc the eigenfunctions of the LO, but not of the 
NLO evolution kernel. The QCD beta- function [3 and the anomalous dimension 7" of the n-th moment a„(p) of the 
wave function have to be taken in NLO [20] . We can substitute the corresponding values into the above equation and 
obtain: 

a 2 (fi b ) = 0.35, a 4 (p 6 ) = 0.18 (LO); a 2 (p 6 ) = 0.218, a 4 (p b ) = 0.084 (NLO), (20) 

at the scale p& = \Jm l B — m' 2 w 2AGeV, which characterizes the mean virtuality of the b quark. The new analysis 
of the experimental data on the j-f*n and it electromagnetic form factor indicates that the twist-2 wave function is 
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close to its asymptotic form [23]. In this article, we use both nonasymptotic and asymptotic form for the w twist-2 
light-cone wave functions and compare the results. 

The subleading twist-4 contributions are presently known only in zeroth order in a s [21,22]. As the twist-3 
contribution is eliminated, we need only the twist-4 wave functions: 

tp ± (ai) = 30S 2 ( ai - a 2 )al[^ + 2e(l - 2a 3 )], <p±(ai) = 30«5 2 ^(1 - o 3 )[^ + 2e(l - 2a 3 )], 

A 



^||(ai) = 1205 e(ai - a 2 )a 1 a 2 a 3 , <p\\ (a*) = -1205 aia 2 a 3 [- + e(l - 3a 3 )], 



/X "2 2 2 

5iN = 2 £ u u 



1 6 

-£(5 2 [uu(2 + 13uu) + 10u 3 lnu(2 - 3m + -u 2 ) + lOu 3 lnu(2 - 3u 

Z 



6 -2M 

5« )], 



92{u) = ^-8 2 uu(u - u). 



(21) 



with S 2 ([i b ) = 0A7GeV 2 and e(/i&) = 0.36. Unlike the case of the twist-2 wave functions, these twist-4 wave functions 
seem to be very difficult to test by experiment, for they usually are of negligible contributions in the sum rules. 

Another important input is the decay constant of B meson f B . To keep consistently, we have to calculate the 
two-point sum rule for f B up to the corrections of order a s . Here we use the two-loop expression for the running 
coupling constant with Nf = 4 and A^ 4 ) = 234 MeV corresponding to a s (Mz) = 0.112 [20] for comparing with 
the results in Ref. [8]. As the value of \x concerned , we take the value 2AGeV which corresponding to the average 
virtuality of the correlation function which is given by the Borel mass parameter M 2 . In the present case a chiral 
current correlator is adopted to delete the contributions from the twist-3 wave functions, we consider the following 
two-point correlator: 



Ti(q 2 )=i J d 4 xe^(0\q(x)(l + l5 )b(x), 6(0)(1 - 75 )<z(0)|0). 
The standard manipulation yields three self-consistent sets of results: 



(22) 



fB 
fB 
fB 



218MeV,m h 
2l2MeV,m b 
206MeV, m b 



4.7GeV, s 
A&GeV, s 
A.QGeV, s 



33Gey 2 ; 
34GeV r2 ; 
35GeV 2 . 



(23) 



The corresponding a s = results: 



fB 
fB 
fB 



l63MeV,m h 
158MeV,m b 
l53MeV,m b 



4.7GeV, so 
A.8GeV, s 
4.9Gey, s 



33Gey 2 ; 
34Gey 2 ; 
35Gey 2 . 



(24) 



From the above results we can see that w 76%, in other word, a, corrections increase the value of fn about 

30%. They will be used as inputs in numerical analysis of the sum rule for fg^iq 2 )- As for the B meson mass mg 
and the pion decay constant f n , we take the present world average value raj = 5.279GeU, and f n = 0.132GeU. The 
continuum subtraction s is about 33 — 35GeV 2 and the pole mass for b quark is taken as m b = 4.7 — A.9GeV. Here we 
make some comments about the continuum subtraction s . The special chiral current leads to cancellations between 
the condensates, the dominating contributions come from the perturbative parts and the nonperturbative parts only 
play tiny roles. The lowest pseudoscalar resonance appears at the energy threshold about s = m 2 B w 28GeU 2 . Though 
the B meson has a narrow decay width, the values taken in Ref. [11] s = 30 — 33GeV 2 are too low due to the large 
difference between the corresponding results for the values of f B . 
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We exploit the sum rule numerically in the following: 







/b/^(0) 


= 60.5MeV,/+ r (0) 


= 


.277, m b 


= 4.7GeV, s 


= 33GeV 2 ; 








= 56.8MeF,/+ w (0) 


= 


.268, m b 


= A.8GeV, s 


= 34GeV 2 ; 






/b/^(0) 


= 53.4MeV, /+,(()) 


= 


.259, m b 


= 4.9Gey,s 


= 35Gey 2 , 


for a s 


^ in LO. 


















/b/^(0) 


= 59.6MeV,/+ T (0) 


= 


.273, m b 


= 4.7Gey,s 


^33GeV 2 , 


for a s 


^ in NLO. 


















/b/^(0) 


= 47.3Mey,/+ w (0) 


= 


.290, m b 


= 4.7Gey,s 


= 33GeV 2 ; 






/b/^(0) 


= 44.1MeV,/+ r (0) 


= 


.279, m b 


= 4.8Gey,s 


= 34GeV 2 ; 






/b/^(0) 


= 41.2MeV,/+ r (0) 


= 


.269, m b 


= 4.9Gey,s 


= 35GeV 2 , 


for a s 


= in LO. 


From the above results we can 


see 


that -7— 

IB 


/+„(0)(a s #0) ^ 
/+ x (0)(a 3 =0) 


s 130%, in 



(25) 



(26) 



(27) 

130%, in other word, a s corrections 
increase the value of /s/g^O) about 30%. Due to the same corrections to the decay constant, the resulting net 
a s corrections are very small, say, for /^.(O) less than 3%. The large correction for /s/g ff (0) is cancelled by the 
corresponding value for Jb- They are compatible with the values obtained in Ref. [8], for a s = 0, fg„(0) = 0.30; for 
a s ^ 0,7^.(0) = 0.27. Our numerical results show that the vibrations for the form factor /^.(O) are about ±0.01 
around the center values, for a s ^ 0, /^(O) = 0.27 ; for a s = 0, /^(O) = 0.28 with LO wave functions. It is shown in 
figure 1 that the form factor fg^q 2 ) with a s corrections lies below the un-corrected one for LO wave function (WF) 
; the quantities of the a s corrections increase with q 2 , at q 2 = 15GeV 2 , numerically lesser than 20% for LO wave 
functions; the curve for NLO wave function lies a little above the corresponding one for LO wave function; the curve 
for asymptotic wave function with a s corrections almost the same as the un-corrected one for LO wave functions at 
q 2 >- 8GeV 2 ; the deviation of the curves for the a s corrected LO wave function and asymptotic wave function from 
each other is notable. In figure 2 and figure 3, we plot the fg^(q 2 ) as function of q 2 in leading order tt light-cone wave 
function with different boundary conditions. From two figures, we can see that the vibrations for /^.(O) are small, 
numerically about ±0.01 around the center values both for the a s corrected and un-corrected form factor. In figure 
4, we use the parameters obtained in Ref. [23] as input, from the figure can see that the curve for f^^iq 2 ) with a s 
corrections varies according to the tt twist-2 light-cone wave functions, the largest deviation of the values from each 
other is less than 15%. In figure 5, we plot the /s T (<? 2 ) with boundary condition s = 33GeF 2 , m b = A.TGeV both for 
a s corrected and un-corrected form factor using the parameters obtained in Ref. [23]. Again, we can see that the net 
a s correction is small. There is a platform for fs^(q 2 ) as function of Borel parameter M 2 for M 2 = 8 — lAGeV 2 which 
verify the value we taken M 2 = 12GeV 2 in calculation. For example, the product fg v (0) is plotted as a function of 
M 2 in figure 6. The uncertainties due to the Borel parameter M 2 can thus be diminished or eliminated. 



V. CONCLUSION 



To summarize, we have re-examined that the weak form factor /^(q 2 ) up to q 2 = WGeV 2 for B decays into 
light pseudoscalar mesons by taking the contributions of a s corrections to twist-2 terms in light-cone QCD sum rule 
framework. Due to the special structure of the chiral current, the contributions of a s corrections to twist-2 terms 
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are of the same Dirac structure as that of Ref. [8] with different weight. As the contributions of twist-3 terms are 
eliminated, the uncertainties due to the twist-3 light-cone wave functions which are not understood as well as the twist- 
2 light-cone wave function are avoided. Furthermore, the possible pollution from wrong parity + mesons are deleted 
by suitable choice of continuum subtraction parameter sq, the final results are supposed to be with less uncertainties. 
The results presented here will be beneficial to the precision extracting of the CKM matrix element \V u b\ from the 
exclusive processes B — > ixivi (I = e,fi), by confronting the theoretical predictions with the experimentally available 
data. Although the a s corrections to SbSb-tx are large, about 30%, the similar corrections to fs canceled them, 
and the resulting net corrections to form factor fg^iq 2 ) are small. Our results are compatible with the observations 
made in Ref. [8] . Compared with the results obtained in Ref. [8] , our results are with lesser uncertainties due to the 
elimination of the twist-3 light-cone wave functions. 
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